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Abstract 

We calculated the energy and the size of the three-deuteron muonic molecule 
(3d2e~ fx") = T>3fj,. It turns out that this system possesses two equilibrium 
positions, one at distances typical for muonic molecules and a second one at 
the usual molecular size. We show, moreover, that the fusion probability of 
the three deuterons is considerably enhanced due to the existence of a 6 Li* 
threshold resonance. Our estimates indicate that this probability is consider- 
ably higher than the decay rate of the competing Auger transition. 
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I. INTRODUCTION 



The physics of two-atomic muonic molecules, as pp/j,, dd/j, or dt/j,, has a long and intense 
history since the prediction of such molecules in the forties and their discovery in the fifties 
Unfortunately, attempts to produce energy for practical use by muon-catalyzed fusion 
in the dtfj, case, the most promising candidate, were up till now unsuccessful. However, the 
corresponding investigations revealed some peculiar aspects of the molecular and nuclear 
physics of such systems. Muonic molecules, e.g., offer the possibility of measuring the 
strong interaction between charged nuclei at very low energies (E ~ eV) which, without 
the Coulomb shielding by fi~, is not accessible to experiment under normal laboratory 
conditions. 

An important aspect of these systems, pointed out in @-^], concerns the formation of thresh- 
old resonances of the nuclei involved, which leads to a considerable enhancement of their 
fusion probability. Such resonances occur in a number of systems, as dtfi, 1 3 He/z or d 6 Li/i. In 
the dtfi case jL| , for example, the near- threshold state 5 He(3/2+) increases the cross section 
for the nuclear transition 

dtfjL — > 5 He(3/2+)/i — > 4 He + /i + n (1) 
by four orders of magnitude as compared to reactions where no threshold resonances occur. 

In the present work we investigate an analogous three-atomic situation, namely, the system 
(3d2e~/i~ = D3/i) in which one of the electrons in the deuterium molecule (3<i3e~) = D3 
is replaced by a muon. Fig. 1. shows that three deuterons can form two 6 Li* resonances 
at energies very close to the 3d threshold energy, Ef£ = 25.32 MeV, i.e., to the energy just 
needed to disintegrate 6 Li into three deuterons. This indicates that in D 3/i three-body fusion 
can take place enhanced by the formation of the (J n , T) = (3 _ , 0) resonance, 

3dfx — ► 6 Li*(3~,0) +fi. (2) 
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It is interesting to note that the D3 molecule and its D3" ion are stable and have quite 
peculiar spectra |7J]. The D3" ion can be produced in the collision 

D 2 + D+ — >D+ + D (3) 

with a high formation probability of the order of A ~ 10 9 cm -3 sec -1 . Thus, it can be easily 
produced under laboratory conditions ||. We also mention that the hydrogen analog of 
the deuterium ion D3", discovered already at the beginning of this century by J.J. Thompson 
P|, plays an important role in the atmosphere of Jupiter and, in general, in astrochemistry 



where it protonates the otherwise chemically inactive CO molecule 10 . 



H+ + CO — > H 2 + HCO+. (4) 

It is, therefore, expected that a sufficient number of D 3m molecules can be produced via the 
D3 production and a subsequent capture, so that there is enough material for the 3d 
fusion process. 

In Sect. II we calculate the energy and the size of the D3/! molecule using a Born- 
Oppenheimer-type anzaz for the wave function. An estimation of the nuclear transition 
rate and, in particular, of the enhancement factor due to the intermediate 6 Li* resonance 
formation is obtained in Sect. III. In order to compare this rate with the magnitude of the 
conventional atomic decay of this molecule, we estimate in Sect. IV the corresponding Auger 
transition probability. Our conclusions are drawn in Sect. V. Finally some technical details 
concerning the Auger transition are presented in the Appendix. 

II. ENERGY AND SIZE OF D 3/ u 

We assume the motion of the heavy particles, the deuterons, to be of adiabatic character 
and therefore their kinetic energy in the total Hamiltonian 

1 IK A X 1 A 1 1 1 

H = -7T- A Pa + A P 2 " + I 1 + I 1 + I 1 

2m e ^ ^ 2 2m M |Pi — P 2 I l r -Pil l r ~ P2I 
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can be omitted. Here g 1 , g 2} and r are the coordinates of the electrons and the muon 
in the center-of-mass (CM) system (see Fig. 2) of the three deuterons. The r-j denote the 
coordinates of the deuterons. 



For the ground state of the D 3 /j molecule we use the anzatz 

^ = X} ^ e -«/2(r - z/ ri ) 2 + e -a/2(r - ur 2 ) 2 + e -a/2{r - z/r 3 ) 2 ^ e -/t/2p 2 - n/2p 2 2 

(6) 

where in the three-deuteron wave function ip(R, x, 9) the following hyperspherical variables 
were chosen instead of the position vectors r^, 

R 2 = x 2 + y\ tan X = -, cosfl = n = (x,9) (7) 

x xy 

where 



Employing the anzatz (^) in a minimization procedure in which the parameters a, v and k 
at fixed values of are varied, we obtain the corresponding energy surface U(R,x,9), the 
effective interaction of the deuterons. Due to its variational origin, this potential will provide 
an upper bound to the exact eigenvalue of H at any pair of angles x an d 9. As can be seen in 
Figs. 3 and 4, the minimum of this function at different values of R is reached for xo — 7r /4 
and 6> = ir/2, i.e., for an equilateral position of the three deuterons. Therefore, the ground 
state for the three-deuterons should appear essentially in this symmetric configuration and 
thus the potential Uq(R) = U(R,xo,9q) can be employed to calculate ip(R,x,9). This is 
supported by the fact that, due to the almost classical motion of the heavy particles, the 
ground state energy will be close to the minimum of the potential surface shown in Fig. 3 



and 4. The dependence on a, v and k on R is shown in Fig. 5 while in Fig. 6 we plot 
Uq(R) versus the hyperradius R. As can be seen, Uq(R) has two minima, the first one at 
"muonic" distances and the second one at atomic distances. The later minimum, obviously, 
corresponds to the case where one of the deuterons in a D 2 molecule, is replaced by the dd/x 
bound system. 

The eigenvalue equation for the 3d states reads 



1 1 d 5 1 A 2 



2m d R 5 dR dR 2m d R 2 uv ') 



+ U (R))il>(R,n)=£ 3 il>(R,n) , (9) 



with 



a2 4 ( 9 . 2 d Id, n d\ 
sm x \uX VX sin oO ad) 

Due to the ^-independence of Uo(R), the solutions of @ can be written as 

1>(r,n) = x{R)y c (n), (11) 

where the hyperspherical harmonics 3^c(^) satisfy the eigenvalue equation 

A 2 y c = c(c + i)y c (12) 

of the grand angular momentum operator A. The quantum number £ = n + Zi + /2 + 3/2is 
composed of the angular momenta l\ and 1% corresponding to the Jacobi coordinates x and 
y, while n denotes the degree of the relevant Jacobi polynomial. 

In a wide range of R (0.1 < R < Rq) and for the quantum number £ = 5/2 (relevant in our 
context) we have 



2m d R 2 

Equation (0) thus reduces to 



<\U (R)\. (13) 



' ' d ; R 5 ^ + U (R))x(R) = £ 3 X(R). (14) 



2m d R 5 OR dR 
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Having in mind the extremely large width of the barrier between the two minima of the 
potential Uq(R) one can consider the states bound by the short-range minimum as physically 
stable with respect to a decay into the second-minimum states. In this area the potential 
can be approximated fairly well by 

U (R) = a + b(R - Rf . (15) 

The constants a and b were found to be (we use the units e = % = m M = 1) 

a = 0.31, 6 = 0.01, i? = 4.45. (16) 

Equation (14) can then be easily solved with the results for the energy and the size of D 3/ u 
being 

£ 3 = -1.74keV, Rq = 11.0 10 -11 cm. (17) 

Both values are of the same order of magnitude as the ones of the corresponding two-atomic 
muonic molecule (£ 2 2.7^ -2.8 keV, Rq ~ 510~ n cm) 0. 

Apart from effects due to the 6 Li resonance formation (see forth section III) comparable 
nuclear transition rates should be expected in both systems. 



III. NUCLEAR TRANSITIONS 

According to Fig. 1 there are two 6 Li* resonances with energies close to the 3d threshold. 
The lowest one for the (J n , I) = (4~, 1) quantum numbers is given by E+iT = (25+i4) MeV 



and the higher one for (J 71 , T) = (3~, 0) [12| which has a broad width E = 26.6 MeV. Since 
the isospin of the deuteron is zero, the 3d state can undergo a fusion reaction only via the 
latter resonance, and hence the D3/X molecule has to be in a J w = 3~ state. 



Following the procedure suggested in |TT| , the probability of the nuclear fusion is obtained 
from the overlap integral M. between the molecular and resonance wave functions, i.e., from 
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w = £ 3 |.M| 2 



(18) 



In contrast to the calculations of the energy and the size of the D 3 /z molecule, where simple 
approximations such as fll^D were sufficient, the calculations of the overlap integral is much 
more delicate. This is due to the fact that the wave function of the three outgoing charged 
particles in the near-threshold region we are interested in, is highly oscillatory Thus, in 
order to avoid direct numerical integrations we choose wave functions that describe the 
crucial input for the D^n molecule and the 6 Li resonance, and allow us to perform the 
integrations analytically. 

For this purpose we approximate the potential Uq(R) at short distances by an analytic 
function. This is easily achieved by using the effective Coulomb potential Z e ff/R. The 
effective charge Z e ff can be found by fitting the eigenpotential U (R) at small distances the 
result being Z e ff = 2.8. In this region the wave function is, therefore, given by the regular 
Coulomb solution Fc(R). 

At large distances the wave function decrease exponentially and therefore we may use a 
phenomenological anzaz of the form 

Vw= ^F £ ( V ,kR)e- kR Y J > M (x,y)u £ (R), (19) 

which describes both the size of the molecule and the Coulomb repulsion of the three 
deuterons at small distances. N m is the normalization constant, rj = Z e ffmd/k is the 
Sommerfeld parameter, and k = \j2md\£z\- Y J ' M (x,y) are the bispherical harmonics cor- 
responding to the total angular momentum J. Since we assume the molecule to be in the 
ground state with negative parity, we have n = and for the sum of the relative angular 
momenta l\ + £2 — 1. In other words we have C = 5/2. 

Due to the repulsion of the deuterons at small distances, the wave function ([19]) is exponen- 
tially small in this region. This means, the main contribution to the overlap integral comes 
from distances R much larger than the nuclear size r$ and thus we may use the asymptotic 



7 



form of the 6 Li* resonance state which is taken to be the asymptotic wave function of three 
charged particles in the continuum 

^res(R, A) = ^f J (R, u)Y J ' M (x, y), (20) 

where 

f J (R,u) = f exp[iqR + ^ ln( 9 fi)]y J ' M (i^) dxdy 
J q 

The factor N res represents the normalization in the nuclear volume and the angles are de- 
noted by Q = (u>, x, y); v(Q) represents the angular part of the total Coulomb potential of the 
three-body systems written in terms of hyperspherical variables while q = ^2m d \E — E th \ 
with E being the energy of the outgoing particle and E t h the threshold energy. 



After some tedious calculations, with r)(Q) = v(Q)m,d/q ~ rj >> 1, and uc(R) = 1, we found 
that the overlap integral, 

M=J ip res (R,Q)ip mol (R,Q) R 5 dRdQ, (21) 



has the behavior 



„ iS(t-) 1 

M ~ e 7T 2 ; (22) 



with 



7 = 7, (J= ^ = 1-06, 

S(t) = [- ln(t) + ln(l -t)]- + ln(l - i( 7 - 1) - 2it) , 

(7 

and 

t = - ^^ + a 2_ 2ia [ 1 + l ], 

fli is obtained from the condition 

Vv(Q) = 

and turns out to be the same as f2o = (xo, $o) introduced above. It should be noted that 
any other reasonable choice of uc(R) is not expected to lead to a significant modification of 
M. 

8 



From (|22|) it follows that an enhancement of the transition rate is achieved if 

ImS(M<0, |ImS(t - )l >l (23) 

7T 2 

In Table | we present the ImS(t_)/7r for certain values of the parameters 7. Evidently the 
condition (^) is satisfied. In other words, the usual decreasing Coulomb barrier is replaced 
by a factor which grows with increasing rj. Since in the present case (near threshold reac- 
tions) rj 3> 1, a considerably enhanced transition of the three deuterons into the 6 Li*(3~, 0) 
resonance is to be expected. The corresponding transition probabilities are given in the 
last column of Table I. The above nuclear fusion process is accompanied by a radiative El 
transition into the 6 Li ground state yielding an energy of about 26 MeV . 



IV. AUGER TRANSITION RATE 

Since the binding energy of D^/j, is essentially higher than the one obtained for the ddfi- 
molecule 0, we have also as a possible decay process the Auger transition which for the 
D3/i molecule is given by 

D 3 „ — > X + e . (24) 

Here X includes the dd\i molecule, the deuteron d, and a second (bound) electron e~. To 
specify the wave function of the system X the following arguments can be used. In the 
Auger transition ( f24"|) most of the energy released is taken by the outgoing electron. This 
means that the relative velocity of the dd/u, molecule and a nucleus d in the final state, is 
much smaller as compare to the velocities of electrons. Furthermore having in mind the 
small distances between all three deuterons in the initial state, one can consider the bound 
electron as moving in the Coulomb field of a "united" atom, formated by the (ddfi) + ion 
and the nucleus d. 

The transition rate for this mode of decay is, of course, quite essential and is of relevance to 
the competing nuclear transition. The corresponding transition probability is determined 
from the amplitude 
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M A =<il>f\V c \il>i>, (25) 

where Vc is the Coulomb interaction of the outgoing electron with the second electron and 
with the center of mass the complex [ddfi + d]. The initial state ijji is the D3/X wave function 
described Sec. II, Eq. (|6|). The final state if)f is chosen to be a product of the wave functions 
of the fragments, namely, 

Ipf = Ipreli'at^out^ddf, (26) 

where VVeZ is the Coulomb wave function for the relative motion of the ddfi molecule and the 
deuteron d, ip a t is the atomic wave function of the bound electron in the field of the ddfi + d 
system, (p out is the plane wave function of outgoing electron, and ipM^ is the ground state 
wave function of the dd/j, molecule which is assumed to be of the form 



ddfi 



? e -Cr d /a^LCAO (2?) 



where i^ddjf designates the Linear Combination of Atomic Atomic Orbitals (LCAO) wave 
function, 

^lcao = 1 1 (e-IW 2 " + e-IW2 + r,\/aA (2g) 

/2(1 + F(r d )) Jna* \ J 



with 

\ a„ 3 \a„J J 

ct M is the Bohr radius for the d\x atom, r d the distance between the deuterons, r M is the muon 
Jacobi coordinate in the ddfi system, and ( = 0.023 is a constant that characterizes the size 
of the ddfi molecule. Due to the comparatively large size of the D3/1 molecule, it is enough 



to use a ddfj, wave function of the form (j27D which provides the correct binding energy and 
has the correct behavior at large distances. 

The width for the Auger transition fl24|) is defined by (see the Appendix) 

T /, D e 4 f^ dx \M A (x,q(x))\ 2 



2 tflQn 3 Jo ~~ q(x 
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where 



x — — , q(x) = 



k " V A 4 

, l -2fi£ 3 

fc = \ T2 ' = 



e 



2/jl£ 3 rn d (2/2 d + m M ) 

K 6m d + m M 

Numerical integration of | A^a| 2 over all final state momenta permitted by the energy- 
momentum conservation yields the decay rate 

A = E w 10 15 sec _1 . (30) 



V. CONCLUSIONS 

Our calculations show that the effective potential energy for the three-deuterons in the D 3 /i 
system possesses two equilibrium positions the first being at a distance characteristic for 
muonic molecules while the second one is of the same order of magnitude as that of the 
usual electronic molecules. The position and value of the "atomic" minimum of the effective 
potential energy imply the following cluster structure of the D 3 /i system at this distance 

*D s »™*D + <**(ddrie> ( 31 ) 

where ^ d is the wave function of the hydrogen atom and ^{dd^)e is a hydrogen-like wave 
function with the ddn + ion in the center of such "atom" . One should bear in mind that such 
a clustering ensures the lowest possible energy of the system D 3 /i. Indeed, as can be seen in 
Fig. 4, this is just the case. It is obvious, that such a state of the D 3 /x system is stable as 
far as electromagnetic transitions are concerned, but it is not stable with respect to nuclear 
fusion in the dd\i subclusters. 

Let's now discuss in brief the possibility of observing the D 3 /i system. One way is to search 
for the Auger electrons with kinetic energy around 1 keV. A second possibility is, of course, 
more exotic and consists in observing the products of the triple fusion and the formation of 
the highly exited state (3~,0) at 26.6 MeV of the 6 Li nucleus. 
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The enhancement factor of the nuclear transition in the D 3 /z molecule, obtained above, 
should be considered as an indication of the possible strong influence of the near-threshold 
nuclear resonance on the molecular system. 
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APPENDIX A: EVALUATION OF THE AUGER TRANSITION 



For the width T the following definition has been used: 

2 J ^ 2kq \ 2p e 2fi D ) 2kq (2vr) 3 (2vr) 3 ' [ ' 



The matrix elements for the Auger transition A4a ( f25|) are given by 
M A = d 3 r D d 3 pd 3 r d d 3 r^d 3 r e ^ ddl ,(r d ,r^ — e 2r V ae 

^ (_? + _!_) ."a' eT» Eel 1 '""?"-"'' (A2) 

where N var is the normalization constant of the variational part of the wave function (^) 
and Td is the Jacobi coordinate of the third deuteron not belonging to the ddp molecule. 
After some analytical derivations, one remains with the following simplified form 

d 3 r d d 3 r D k 3/2 e~^ r<i 



KA fu . 162^VC /a,y/ 2 r 



/3(l + 2e X p(-i^)) a3/4 y/^ + FM) 



where 



E^>)^ (-^/B)g)^ /K , + x( K/2 ,^)) *(«) (A3) 

ci 2 

1 o 2 o ™ 4tt / 1 



2 3 Varfii VV" 



2n 
4 



- a) e"« - -)V2 erfc f ^ - + a) e"« + er fc f ^±5 



v/2 



v/2 



00 f_^n-l 



3 2 V ; 61 (2n-l)! 



/"oo _ 2/1 

T(a,5) = / dx 2 sinf&r) e aa; ( — JAx, a) + hix, a 
Jo \x 
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Ji(x, a) = e 



1/a f 



l2a 3 / 2 



7(3/2, a{x + 1/a) 2 ) - 7 (3/2, 1/a) 



-1/a _ p —a[x + 1/a)' 



1 /F 



2a 2 V a 



1 — erf (y/a(x + 1/a] 



JJ X) a ) = eV« -1 f e-«( x + V«) 2 _ /IT h _ er f f yfa(x + 1/a) 
2a \ V a L v 



and 



We note here that in ( |^3| ) the following dimensionless variables aa 2 — > a, Ka 2 , 
qa^ -> g, r d /a M -> r d , and r D /a M -> r D are used. 
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TABLES 



TABLE I. The function ImS(t_) for three different values of of the parameter 7 and a = 1.06. 

7 r)(Clo) ImS(t_)/7r w s -1 

0.857 21.20 -0.5244 5.42 10 19 

1.916 9.48 -1.5130 3.37 10 44 

2.710 6.70 -2.2590 3.03 10 50 
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FIGURES 
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6 Li 



FIG. 1. The 6 Li nucleus spectrum. 
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FIG. 2. Jacobi coordinates for all constituent particles of the D3/U molecule. 
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0.70 



FIG. 3. Two-dimensional cross section of U(R,x,&) with R = 4.5 a^. The surface has a 
minimum at 6 = ir/2 and x = 7r /4 which corresponds to the equilateral triangle formed by the 3d 
nuclei. 
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FIG. 4. Same as in Fig. 3 with R = 2.5 a 
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R/a, 



FIG. 5. The dependence of the parameters k, a, and v of Eq. (6) on R 
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FIG. 6. The Uq(R) = U (R, xo, #o) effective potential. On the upper figure the potential is shown 
at atomic distances while on the lower at mesonic distances. The a e and a M are the electronic and 
muonic Bohr radii respectively. 
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